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This paper presents a novel energy management strategy (EMS) to control an isolated microgrid powered 
by a photovoltaic array and a wind turbine and equipped with two different energy storage systems: 
electric batteries and a hydrogen production and storage system. In particular, an optimal scheduling 
of storage devices is carried out to maximize the benefits of available renewable resources by operating 
the photovoltaic systems and the wind turbine at their maximum power points and by minimizing the 
overall utilization costs. 

Unlike conventional EMS based on the state-of-charge (SOC) of batteries, the proposed EMS takes into 
account the uncertainty due to the intermittent nature of renewable resources and electricity demand. In 
particular, the uncertainties are evaluated with a stochastic approach through the construction of 
different scenarios with corresponding probabilities. The EMS is defined by minimizing the utilization 
costs of the energy storage equipment. 

The weather conditions recorded in four different weeks between April and December are used as case 
studies to test the proposed EMS and the results obtained are compared with a conventional EMS based 
on the state-of-charge of batteries. The results show a reduction of utilization costs of about 15% in 
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comparison to conventional SOC-based EMS and an increase of the average energy storage efficiency. 


© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 


According to the IEA [1], today more than one billion people 
worldwide lack access to electricity. Various options for supplying 
this electricity need to be considered: these include on-grid, mini- 
grid, and isolated off-grid configurations. In particular, decentral- 
ized options are an important alternative in cases where grid 
extensions are too expensive. Stand-alone power systems using 
fossil fuels suffer from the high costs of these fuels, pollution and 
greenhouse gas emissions, while stand-alone microgrids com- 
pletely powered by renewable energy sources (RES) can overcome 
these problems. One of the main drawbacks of using energy con- 
version systems based on RES, such as photovoltaic systems (PV) 
and wind turbines (WT), is their discontinuity in energy produc- 
tion, so that an effective energy storage section is required to 
ensure a continuous electrical energy supply. 

Among energy storage systems, batteries are the most common 
choice for short-term storage. However, for longer-term storage of 
energy, their application might be inappropriate owing to their low 
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energy storage density and unavoidable self-discharge [2]. Energy 
storage systems based on hydrogen technologies are one of the 
most interesting options and have been the subject of many stud- 
ies and research activities in recent years [2-14]. In hydrogen stor- 
age systems, excess electricity can be converted to hydrogen 
through an electrolyzer (EL) and stored in pressurized tanks. The 
stored hydrogen can later be used to produce electricity in a fuel 
cell (FC). 

Hydrogen technologies are not yet competitive with large scale 
energy storage systems, such as pumped hydro storage. However, 
they have major advantages in stand-alone power systems. Owing 
to the absence of the electrical grid, stand-alone systems usually 
rely on internal combustion engines (ICE), which use fossil fuels. 
The substitution of these fuels with RES, mainly wind turbines 
and photovoltaic modules, leads to significant reduction of pollut- 
ing emissions. In this framework, the coupling between hydrogen 
storage systems and electrical batteries is investigated as a suitable 
solution in stand-alone power generation plants, which may lead 
to greater flexibility in following load demand and ensuring both 
long and short-term energy storage. Depending on the situation, 
the surplus of electrical energy can be used to charge the electric 
batteries or the hydrogen tanks. Similarly, the demand can be 
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Nomenclature 

Symbols 

p probability 

A panel area (m?) 

C cost (€) 

e percent error 

f objective function 

G solar irradiation (W/m?) 
L lifetime (h) 

N number of elements 
n molar flow (mol/s) 

P power (W) 

p pressure (Pa) 

Q capacity (Ah) 

S number of scenarios 
T temperature (K) 

t time (h) 

U voltage (V) 

V volume (m°) 

v wind speed (m/s) 

Z height (m) 

Zo roughness height (m) 
a temperature coefficient 
n efficiency 

T branching time (h) 

0 time horizon (h) 
Subscripts 

AMB ambient 

B battery 


covered by using the energy stored in one of these two systems. 
Such highly integrated energy systems require a proper energy 
management strategy (EMS) to ensure optimal operation of the 
stand-alone power plant. 

The common control parameter used in conventional EMSs is 
the battery State-of-charge (SOC) [8-11]. Generally, these strate- 
gies do not consider the utilization costs and lifetimes of the differ- 
ent devices. A more accurate techno-economic optimization of the 
overall system requires the inclusion of investment costs, operat- 
ing costs and lifetimes of the main plant components. In particular, 
Dufo-Lopez et al. [12] proposed a control strategy based on the 
minimization of operational costs of stand-alone hybrid systems 
with hydrogen storage. A different EMS, based on a stand-alone 
power system that takes into account the utilization cost of the 
energy storage system, was introduced by Torreglosa et al. [13]. 
The comparative study carried out by Castañeda et al. [14] high- 
lighted that the EMS based on the minimization of utilization costs 
leads to a higher utilization of the hydrogen system compared to 
EMS based on the use of the battery SOC. 

The optimal scheduling of the different microgrid components 
is widely studied in literature. In particular, generation scheduling 
can be divided into two different problems: Unit Commitment (UC) 
and Economic Dispatch (ED) [15]. The UC strategy defines the on/ 
off status of generators over a daily or weekly time horizon while 
respecting system constraints. The ED strategy defines the operat- 
ing power of the units committed by the UC problem for a shorter 
time horizon (hourly or in real-time). The determination of the 
optimal EMS greatly benefits from the inclusion of wind and solar 
radiation forecasts. In particular, Logenthiran and Srinivasan [16] 
proposed a three-step algorithm to solve the UC problem for a 
stand-alone microgrid based on RES and batteries. Zein Alabedin 


D demand 

EL electrolyzer 

EX excess power 

FC fuel cell 

G solar irradiation 
H2 hydrogen 

LD load 

MIN minimum 

MAX maximum 

NOM nominal 

PV photovoltaic panel 
REF reference 

S scenario 

UN undelivered power 
w wind 

WT wind turbine 
Acronyms 


SOC state-of-charge 

NOCT nominal operating cell temperature 
PEM Proton Exchange Membrane 

EMS energy management strategy 

RES renewable energy sources 

PDF probability distribution function 
O&M Operation and Management 

LHV lower heating value 

MPPT Maximum Power Point Tracking 
GAMS General Algebraic Modeling System 


et al. [17] introduced a different methodology based on a multi- 
scenario stochastic model and the assignment of an additional 
spinning reserve requirement to take into account forecast uncer- 
tainties. In particular, two optimization models to solve the UC 
problem in a microgrid including RES and batteries operating in 
grid-connected mode and isolated mode were presented in [17] 
and several cases studies were examined to investigate the opera- 
tion of a microgrid. Tuoy et al. [18] examined the benefits of using 
a stochastic approach instead of a deterministic one to account for 
the uncertainty of wind in the UC problem for a power generation 
system with large wind penetration. Possible future wind and load 
profiles were represented with a scenario tree. The optimal sched- 
uling of a renewable stand-alone microgrid that includes a fuel cell 
was introduced by Morais et al. [19]. The analysis was carried out 
with a deterministic approach for a period of one day and the fuel 
cell was used only as an emergency reserve. A similar analysis was 
presented by Khodr et al. [20], who simulated a renewable micro- 
grid to propose a method for the optimal management of the next 
week (672 time steps) in a deterministic environment. 

All of these works in the field of EMS are based on a stochastic 
approach applied to stand-alone microgrids without hydrogen 
storage systems. This paper proposes a novel EMS that optimizes 
the scheduling of RES stand-alone power systems that include 
batteries and hydrogen storage technologies. To demonstrate its 
capabilities, the proposed EMS is applied to a microgrid powered 
only by RES and currently under construction at the Concentrating 
Solar and Hydrogen from Renewable Energy Sources Laboratory as 
part of the Sardegna Ricerche Renewable Energy Cluster. In partic- 
ular, a mathematical model to minimize the utilization cost is elab- 
orated. A stochastic approach instead of a deterministic one is 
chosen and the scenario tree method is used to take into account 
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the uncertainties of both forecast weather and load profiles. In this 
paper the proposed EMS is presented and compared with a 
conventional SOC-based EMS with reference to a stand-alone 
micro-grid powered by a photovoltaic array and a wind turbine 
integrated with electric batteries and a hydrogen storage system. 


2. Modeling of the hybrid system 


The purpose of this study is to develop an optimization model 
for load generation scheduling of a hybrid micro-grid taking into 
account the uncertainty of both resource availability and energy 
demand. The focus is on a stand-alone hybrid microgrid com- 
pletely powered by photovoltaic systems and wind turbines. In 
particular, the system configuration considered here refers to the 
microgrid of the Hydrogen from Renewable Energy Sources 
Laboratory. The purpose of this laboratory is the implementation, 
testing and demonstration of technologies related to the produc- 
tion, storage and use of hydrogen from renewable energy sources. 
Currently, laboratory activities include the development of proto- 
type plants, the characterization and integration of individual com- 
ponents and the overall power generation system, as well as the 
identification of optimum design criteria in accordance with com- 
ponent specifications and the electricity demand curve. For this 
reason, the laboratory includes a microgrid powered exclusively 
by RES to satisfy the energy demand of the laboratory itself. 
Fig. 1 shows the configuration of the stand-alone hybrid system. 
In particular, the hydrogen storage system includes a PEM electro- 
lyzer, four independent tanks having an overall storage capacity of 
55 N m? (13.8 bar and 1m? each) and a 5 kW PEM fuel cell. To 
study and compare different storage options, a bank of 32 lead-acid 
batteries is also included. Table 1 shows the main characteristics of 
the microgrid devices whereas Table 2 shows their costs and 
expected lifetimes. 

To enable the EMS to predict the performance of the system, a 
simple simulation model of each device is introduced together 
with the assessment of respective utilization costs. In the following 
subchapters, the simulation models for each component are 
described in detail. Despite the lower accuracy, the adoption of 
simple simulation models is essential to reduce the computational 
time and to allow a fast running of the EMS during the operating 
phase of the microgrid. On the contrary, more accurate simulation 
models should be used for the design of the power generation 


system as well as during the verification phase of the EMS 
performance. For this reason, the simulation models used here 
for evaluating the expected performance of the microgrid as a 
result of the generation scheduling defined by the EMS. The latter 
simulation models were already presented in a previous work [26] 
and just their basic features are reported in this paper. 


2.1. Photovoltaic system 


It is assumed that the PV system always operates at its 
maximum power point and therefore its operation is independent 
of the EMS. At a given irradiation condition and module tempera- 
ture, the working voltage of the PV is adjusted by the MPPT (Max- 
imum Power Point Tracking) control system to maximize the 
power output of the array. The following relation is used to 
describe the power output Ppy of the PV system: 


Ppy = G- Arv - Nov - Npy (1) 


where G is global sun irradiance (W/m7), Apy the panel area (m7), 
Npy the total number of panels. The efficiency of the solar panel 
pv is expressed as a function of irradiance G and ambient 
temperature Taye as follows: 

Tre) (2) 


where ‘Ypyvrer is the efficiency of the panel under reference 
conditions (1000 W/m? of solar irradiation and 25 °C of cell temper- 
ature), « is the temperature coefficient (1/K), NOCT is the nominal 
operating cell temperature and Tper is the reference module temper- 
ature (25 °C). All of these parameters are provided by the panel 
manufacturer [21]. 


NOCT — 20 
"py = "pv REF i X(T +G 300 


2.2. Wind turbine 


Similar to the PV system, the power produced by the wind tur- 
bine is independent of the EMS. The power curve is provided by the 
manufacturer and is depicted in Fig. 2 as a function of wind speed 
at hub height [22]. Wind speed is measured at the anemometer 
height (Zan) and therefore a wind shear is modeled to calculate 
the wind speed at the turbine hub height (Zwr). Wind shear is a 
function of both altitude and roughness of the ground. The relation 


Fig. 1. Configuration of the stand-alone hybrid systems. 
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Table 1 
Main components and characteristics of the microgrid. 


Photovoltaic system 


Hydrogen storage system 


Panel peak power 0.225 kw Electrolyzer 

Nominal efficiency, "pv,rer 18.1% Nominal power 6 kw 

Number of panels, Npy 36 Hə production rate: 1.05 Nm?/h 

Panel area, Apy 1.244 m? Delivery pressure 13.8 bar 

Solar cell material Monocrystalline silicon Electrolyte PEM 

Temperature coefficient, « Fuel cell 

NOCT —0.38%/K Nominal power 5 kw 

45°C Hə rated consumption 3.90 Nm3/h 

Electrolyte PEM 

Wind turbine Batteries 

Rated power @ 14 m/s 3 kW Nominal voltage, Ug 12V 

Configuration 3 Blades, vertical axis Rated capacity (C20), Qs 240 Ah 

Rotor diameter 4 Number of batteries, Ng 32 

Hub height, Zwr 5.8 m Batteries per string 4 

Anemometer height, Zan 2m Type Lead-acid 

Table 2 1 T T ii T T T 
Costs and estimated lifetime for storage devices. : 
Batteries Electrolyzer Fuel cell ; 

Investment cost 400 €/battery 75,000 € 28,000 € oa SR 

Estimated lifetime 1300 cycles (DOD = 30%) 30,000 h 30,000 h : A N 

O&M cost - 0.2 €/h 0.2 €/h 


Pyr KW] 


Wind speed [m/s] 


Fig. 2. Wind turbine power curve (Ropatec Simply [17]). 


used to calculate wind speed at hub height (vwr) as a function of 
wind speed measured by the anemometer (van) is: 


Vvwr = Van: In (22) /im =) (3) 


Zo being the surface roughness assumed here as 0.0025 m. Based on 
wind speed at hub height, the power output of the wind turbine is 
calculated from the manufacturer’s power curve. 


2.3. Batteries 


The use of batteries to compensate for the differences between 
energy production and energy demand is the most common 
solution in microgrid projects involving RES. The addition of a 
hydrogen storage system leads to a different managing of batteries 
to maximize their efficiency and lifetime. The available energy 
content of an electric battery is commonly measured through the 
so-called “state of charge” (SOC). The latter is the ratio between 


Th ps (SOC=60)~ ~ 


Efficiency 


—B— Nig cy (SOC=80) 
—2— ngns (SOC=80) 
ty 


‘i mall cc 
o 1000 2000 3000 4000 5000 6000 
Electrical Power [VV] 


Fig. 3. Batteries, electrolyzer and fuel cell efficiency as a function of power and SOC 
(the efficiency of batteries varies within the shaded area according the SOC). 


the stored energy and its nominal storage capacity and is calcu- 
lated by monitoring the charging (Pg,cu) and discharging power 
(Pg.pis) over time: 


(Paci -NB cH T Pais/Naps) At 
NpUsQz 


where "ecu and Ng pis are the battery efficiencies during charge and 
discharge processes respectively, At is the applied time step (h), Ng 
is the number of batteries, Qg is the battery nominal capacity (Ah) 
and Us is the battery nominal voltage (V). Battery efficiency is 
expressed in function of the SOC and the input/output power and 
is calculated in accordance with [23]. The inverter/rectifier has 
the task of controlling and managing the battery bank. The manu- 
facturer provided the inverter efficiency curve as a function of 
power [24]. The battery efficiency including inverter losses is shown 
in Fig. 3 as a function of power and SOC (shaded area). 


SOC(t) = SOC(t — 1) 4 (4) 


2.4. Hydrogen storage system 


The use of hydrogen storage systems, alone or in addition to 
batteries, has been the subject of several studies and research 
activities in recent years. In hydrogen storage systems, an electro- 
lyzer produces hydrogen when the power delivered by RES is 
higher than demand. Applying Faraday’s law, the hydrogen molar 
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flow of the electrolyzer (ny, ¢,) can be expressed as a function of the 
supplied electric power (Pg): 


Ne Pew 
Np, EL = LHV, (5) 


where LHV, is the lower heating value of hydrogen (240 MJ/kmol) 
and ne is electrolyzer efficiency, which takes into account electro- 
chemical, thermodynamic and ancillary losses. 

The hydrogen can be used to produce electrical energy by 
means of a fuel cell during peaking periods. Hydrogen consump- 
tion of the fuel cell (ny, rc) is directly related to its power output 
(Prc) through the following relationship: 


_ Pr 
Np, FC = eclHVi, (6) 


where rc is fuel cell efficiency, which takes into account electro- 
chemical, thermodynamic and ancillary losses. In accordance with 
[25,26], both ye, and nrc are expressed in function of input/output 
power, as shown in Fig. 3. 

An important control variable of a hydrogen storage system is 
the hydrogen tank pressure pp2(t) at a certain time t. Tank pressure 
is a measure of the hydrogen content of the vessels and can be 
expressed in function of the produced (np2 z1) and consumed 
hydrogen molar flows (np2 rc) as well as the pressure at the 
previous time step pp2(t — 1) as follows: 


RT 
Pu, (t — 1) 4 v (Np, EL — Np, Fc) (7) 
H2 


Pu, (t) 


where R is the gas constant, Tp, is the mean temperature inside the 
vessel (assumed here constant and equal to 313 K) and Vy, is the 
overall tank volume (4 m? in this study). It is worth noting that 
the maximum hydrogen pressure is relatively low (13.8 bar) and 
this is why the ideal gas law is used. 


2.5. Excess and undelivered power 


Load shedding and power curtailment are often required to 
maintain the power balance of the system. If demand exceeds 
the capacity of renewable generators and the storage system is 
unable to cover this deficit, load shedding is needed to preserve 
the power balance. Vice versa, if the entire demand can be supplied 
by the renewable generators and storage systems are fully charged, 
any excess of renewable energy cannot be stored and a reduction 
of the power produced by the generators is necessary. Both load 
shedding and power curtailment options are introduced in the 
power balance through two variables: undelivered power (Pyn) 
and excess power (Pex) respectively. 


3. Problem formulation 


Since the system is solely powered by non-controllable renew- 
able sources, PV and WT are operated at all times at their 
maximum power point. Therefore, power output is directly a func- 
tion of solar radiation and wind speed respectively. In particular, 
starting from the forecast of solar radiation, ambient temperature 
and wind speed conditions, the expected power produced by the 
photovoltaic system and wind turbine were calculated for each 
time step At of the overall time period 0. For given values of the 
forecasted load demand, the purpose of the EMS is to determine 
the optimal operation of storage devices (batteries, electrolyzer, 
fuel cell) to give the minimum operating cost of the overall micro- 
grid. According to the scheduling determined by the EMS, the 
power values of each device (P;) are linked with binary variables 
(Yi) that determine the status of the device, as better explained 
in the following. 


3.1. Utilization costs 


When excess power produced by PV and WT generators needs 
to be stored, the EMS allows assessment of the most convenient 
storage system. In particular, the EMS calculates and compares 
the utilization costs of batteries and the hydrogen storage system 
for cycling the energy corresponding to available power P during 
a time step At. Because of the absence of fuel costs, utilization cost 
C takes into account only depreciation and replacement costs Cry 
and operating and maintenance cost Cog, for each device used: 


1 Ci i 
=a ( L + Coumi) (8) 


where L; is the lifetime of the i-th device and y is the roundtrip 
efficiency. The latter has been included in Eq. (8) to take into 
account the different energy losses inside the storage systems and 
to carry out the comparison referring to the same output and not 
input energy. 

Starting from Eq. (8) and according to [12], the battery utiliza- 
tion costs Cg, cą during charging are calculated as: 


Cegn /Le,cu + Coem B 
"Ne.cu"ls pis 


CBH 


(9) 


where Cp is the capital cost of the battery bank (€) and Lg cy (h) is 
the battery lifetime evaluated during charging while the product of 
the battery efficiencies during charging and discharging processes 
gives the battery roundtrip efficiency. The batteries’ O&M costs 
Coam,.s are of minor importance and are neglected in this study. 

The lifetime Lgcy in terms of hours is difficult to know, while it 
is more significant to refer to the equivalent full charge/discharge 
cycles Neycies 


Neyctes (1 0) 


In Eq. (10), the ratio of the cycled energy and its nominal capac- 
ity gives the battery operation time in terms of equivalent full 
cycles. 

Similar to battery utilization costs during charging, the hydro- 
gen storage utilization costs refer to the costs of producing hydro- 
gen through the electrolyzer and using it for fuelling the fuel cell: 


(Cerin /Ler + Coemt) + (Crcwn/Lec + Cosme) : Atrc/At 
Ner ` Nrc 


Ch2,cH 
(11) 


where Cez in and Crcyn are electrolyzer and fuel cell acquisition costs, 
Le and Lec electrolyzer and the fuel cell lifetimes, Coguer and Cogmrc 
the O&M costs of the electrolyzer and the fuel cell, ng, and nrc are 
the electrolyzer and fuel cell efficiencies and Atpc/At is the opera- 
tion time of the fuel cell. Unlike the assumption of [12], the fuel cell 
operation time is different from the electrolyzer operation time. The 
operation time Atrc/At is calculated assuming that the fuel cell 
operates at its nominal power Pomc and uses all the hydrogen pro- 
duced by the electrolyzer during the time step At: 


Ate /At = Ne, re p= (12) 
'NOM,FC 

Similarly, if the user’s energy demand has to be provided by the 
storage system, the EMS calculates the cost of supplying the energy 
with the batteries or the fuel cell. 

If the energy demand has to be provided by the storage system, 
battery utilization costs during the discharge phase are equal to 
the average cost of supplying a certain power Pz pis with the batter- 
ies for a time step At and is calculated as: 


Cz pis = Coin /Le,pis a Coxm.s (13) 
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where Lg ps (h) is the battery lifetime evaluated in a discharging 
step and defined as: 


NgUBQg 
Le pis = 5——— Nevcies (14) 
Pe.pis/Np.pis 
If the energy demand has to be provided by the storage system, 
the utilization cost of supplying the energy with the fuel cell Cgc is 
the average cost of supplying the power Prc for a time step At: 


Cre = Cren /Lre + Cogeme (15) 


Finally, to avoid load shedding and power curtailment as much 
as possible, virtual costs were associated with undelivered power 
(Cun) and excess power (Cex). In particular, the specific costs for 
undelivered and excess power were set to 0.005 €/W h. 


3.2. Scenario tree 


The scenario tree approach was adopted to introduce uncertain- 
ties into the methodology, since weather data can be forecast only 
within a very limited time frame and is always subject to devia- 
tions. Similarly, another system variable that cannot be directly 
controlled is the actual energy demand of the microgrid, which 
depends on many factors (weather conditions, occupants’ behav- 
ior, etc.). To deal with these uncertainties in energy consumption 
and electricity production from PV and WT, a stochastic approach 
was chosen. Tuohy et al. [18] compared stochastic and determinis- 
tic methods for scheduling energy systems with high wind pene- 
trations and state that stochastic optimization models perform 
better in terms of operational costs and precision. In this work, a 
scenario analysis is applied as a special stochastic method. 

Forecasting errors in load profiles, wind speed and solar radia- 
tion are usually represented through a probability distribution 
function (PDF) curve as shown in Fig. 4. The generation of the dif- 
ferent scenarios is carried out through a discretization of the PDF 
curve of the forecasting errors by a set of finite states (solid dots 
in Fig. 4) such that a probability is assigned to each state [17,28]. 
In this way, by coupling the error value of the ith state (e')with 
its corresponding probability (p’), it is possible to define discrete 
probability distribution sets for load demand (6p), wind speed 
(ôw) and solar radiation (ôç) in accordance with [17]: 


dp = {(€p, Pp); (€D; PB): «+++ (€D Po) } (16) 
dw = {(Cw, Pw): (Ew Pw): (ws Pw) } (17) 
dc = {(€6, Pa), (CG pè), -- (€8, pe), } (18) 


where n, m and q are the number of states in the discrete set of load, 
wind speed and solar radiation forecasting errors. A set of scenarios 
can be created from the discrete sets dp, dw and ôç to take into 
account possible deviations from the load, wind speed and solar 
radiation forecasting values. The product of the number of states 
n, m and q is equal to the total number of scenarios (S), whereas 
the probability for each scenario (p) is equal to the product of the 
probabilities of the states corresponding to that scenario. 

As described in [27], the usual approach for a scenario analysis 
is to model the set of S independent scenarios to find the best one. 
However, this approach is very calculation intensive, since in many 
cases a large number of scenarios have to be considered. Moreover, 
the calculations must be continuously replicated when the new 
forecasts are available. To reduce the number of decision variables 
Pallottino et al. [27] used a so-called scenario tree approach, which 
aggregates into a single bundle the set of scenarios sharing a com- 
mon portion of data. In this way, the common portions of the S 
independent scenarios are aggregated into bundles by producing 
a tree structure, as shown in Fig. 5. The main concepts behind 
the aggregation process of the scenario tree are the branching time 
and the stage. The branching time t is the time in which the 
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Fig. 4. Discretized probability distribution functions of the load, wind speed, and 
solar radiation forecasting errors. 


scenarios begin to differ, while a stage is the time period between 
two branching times. Each dot of Fig. 5 represents the state of the 
microgrid during a time step At. Up to the first branching time, 
that is during stage 0, all decisions for the different scenarios are 
the same. After the first branching time, stage 1 begins with differ- 
ent decisions for each scenario, and so on through the following 
stages. During stage 0 of the scenario tree the state of the microgrid 
is common to all the scenarios. The main rules adopted to organize 
the scenario tree are described in [27]. 

Rolling planning is used to update data and permit a continuous 
revision of the results. Starting from the current conditions, the 
storage systems are scheduled until the first branching time. This 
schedule is updated for any given rolling time. 


3.3. Mathematical model 


The purpose of the EMS based on short-term generation sched- 
uling is to determine for each time step the status of the overall 
microgrid giving the minimum operating costs. For this reason, 
the objective function (f) is formulated in such a way as to 
minimize the sum of the utilization cost of devices (batteries, 
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Fig. 5. Scenario tree structure. 
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electrolyzer and fuel cell) included in the two energy storage sys- 
tems during the charging/discharging processes and costs related 
to excess and undelivered power: 


0 


s 
f = min X` p(s)X _[Cscn(s, t) + Ce pis(S, t) + Cuzin (S, £) + Crc(S, t) 
s=1 


t=1 


+ Cex(S, t) + Cun(S, t)] (19) 


The minimization problem is subject to the following 
constraints: 


(1) Energy balance equation: 


Ppy(t) + Pwr(t) + Ps pis(t) - Ys pis + Prc(t) - Yre + Pun(t) 
= Pip (t) + Pacu(t) - Ys cu + Pe(t) - Yer + Pex(t) (20) 


where binary variables (Y;) determine the status of the device 
(1 =on, 0= off). The overall input power in the microgrid has to 
be equal to the overall output power, whatever the time step and 
scenario. 

(2) Unit generation limits: 


Pgcu(t), Pa.pis(t) < Pamax (21) 
Permin < Pelt) < Permax (22) 
Prcmin < Prc(t) < Promax (23) 


For batteries, maximum charging and discharging power 
(Ps.max) Was introduced since efficiency decreases considerably 
for high power flows. For the same reason, a limit on maximum 
power was introduced for electrolyzer (Permax) and fuel cell 
(Pecmax). In addition, although the electrolytic cells can efficiently 
operate even at low power levels, a significant increase in ancillary 
losses occurs at low values of power with a decrease in efficiency. 
Therefore, a limit on the minimum power of electrolyzer and fuel 
cell (Per mın and Prcmin) Was also imposed. 


(3) Storage limits: 


SOCuin 


PH, MIN 


SOC(t) < SOCmax (24) 
Puo(t) < Puh, max (25) 


The batteries usually have a minimum SOC (SOCmm) 
recommended by the manufacturer below which they should not 
operate. The minimum SOC value greatly affects battery lifetime 
because a deep discharge of the batteries causes a significant 
decrease in the number of cycles before replacement. The maxi- 
mum SOC (SOCmax) is usually reached when batteries are fully 
charged (100% of the SOC). 

Since a compressor is not included in the hydrogen storage 
section, the maximum pressure of the hydrogen tanks (py, max) 
corresponds to the hydrogen delivery pressure of the electrolyzer. 
In addition, a constraint on fuel cell supply pressure is imposed by 
the manufacturer. This limit therefore corresponds to the mini- 
mum pressure achievable in the hydrogen tanks (Py, mn). 


< 
< 


(4) Congruity constraint: 


Yperpc(t < T,S1) = Yperrc(t < 7,52) VS1,S2 ES (26) 


This constraint requires that the binary variables determining 
the status of batteries, electrolyzer and fuel cell in each scenario 
of the same bundle are identical up to the following branching 
time. 

The formulated problem involves continuous and binary 
variables. The nonlinear efficiency curves of batteries, electrolyzer 
and fuel cells introduce a nonlinear behavior of utilization costs 
(Eq. (9), (11) and (13)) as well as of SOC (Eq. (4)) and pu (Eq. 
(7)). However, through a piecewise linear approximation [29] it 


is possible to solve the problem as a mixed-integer linear program- 
ming (MIP) problem using Gurobi [30]. 


4. Case study 
4.1. Input data 


The capabilities of the optimization approach presented were 
evaluated with reference to the microgrid configuration of the 
aforementioned Hydrogen from Renewable Energy Sources Labo- 
ratory. The analysis was carried out with reference to four different 
case studies: a week in August (summer case), a week in October 
(autumn case), a week in December (winter case) and a week in 
April (spring case). Time step At and branching time t were set 
equal to 0.25 h and 6h respectively and a time horizon of 4 days 
was set to take into account the possible evolution of load and 
weather. All these parameters were subjected to a sensitivity anal- 
ysis and were optimized for the physical problem under consider- 
ation. Forecast data of solar radiation, wind speed and direction as 
well as ambient temperature with a sampling time of 1 h were pro- 
vided by accredited institutions [31,32] and refer to the weather 
station in the ENEA Casaccia Research Centre (near Rome). A roll- 
ing time of one hour was imposed and weather forecast data were 
updated every 24h. Starting from the forecast of solar radiation 
and wind speed conditions, the expected power produced by the 
photovoltaic system and wind turbine were calculated. This 
expected power represents the main scenario, but other scenarios 
were introduced with a corresponding probability to take forecast- 
ing errors into account. Concerning the load profile, since mea- 
sured data were not available, according to the information 
available on laboratory activities, a variable load with an 8 h peak 
period centered at noon was assumed. In particular, a peak load of 
3200 W and a base load of 800 W for the remaining hours were 
assumed with an overall daily energy demand of 38.4 kWh. As 
demonstrated in a previous paper [34], the latter load values max- 
imize the energy supplied to the microgrid and ensure that the 
demand is covered throughout the year. 

Data for the PDF of forecasting errors in load profiles, wind 
speed and solar radiation were extrapolated from [28]. The PDFs 
for the load profile and wind speed were discretized into 5 states 
whereas only 3 states were used to discretize the PDF curve of solar 
radiation, as shown in Fig. 4. Consequently, a total of 75 scenarios 
were obtained. To decrease computational complexity and running 
time, the total number of scenarios was reduced to 15 by applying 
a reduction algorithm as described in [33]. The difference between 
the most extreme scenarios in terms of production and demand of 
energy is about 30%. 

Table 3 shows the values imposed for the initial conditions and 
maximum and minimum storage limits. 


Table 3 

Limited and initial value imposed. 
Battery 
Minimum state-of-charge (SOCyin) 60% 
Maximum state-of-charge (SOCyax) 90% 
Initial state-of-charge (SOC) 80% 
Maximum charging/discharging power (Ps,max) 18 kw 
H2-storage 
Minimum H; tank pressure (py wn) 2 bar 
Maximum H; tank pressure (Dy max) 13.8 bar 
Initial H3 tank pressure (py Jv) 10 bar 
Minimum electrolyzer power (Pg min) 1.5 kW 
Maximum electrolyzer power (Pg max) 6.2 kw 
Minimum fuel cell power (Prc.muv) 0.5 kw 
Maximum fuel cell power (Pfc max) 6 kW 
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The current case study involves 18367 variables of which 12360 
binary variables for each scenario tree built at a certain rolling 
time. 


4.2. Results 


Starting from the initial conditions, the solver gives the solution 
of the UC problem, which is the optimal scheduling of the various 
devices that minimizes total utilization costs. The main output of 
the program is therefore the on/off status of the various devices. 
Fig. 6a-d shows the optimal schedules for the electrolyzer and fuel 
cell for the four weeks studied. As shown in Fig. 6a, the electrolyzer 
is widely used during the summer week. In this period, a surplus of 
energy occurs during the central hours of the day mainly due to high 
PV energy production. The deficit of energy at night is instead sup- 
plied by the batteries while the use of the fuel cell is not required. 

A less frequent use of the hydrogen storage system occurs dur- 
ing the autumn week (Fig. 6b) and the spring week (Fig. 6d). In 
these cases, the batteries are mainly used to compensate for the 
minor difference between energy production and demand, whilst 
the hydrogen storage system works only for few hours. In contrast, 
Fig. 6c shows a wide use of the fuel cell during the winter week, 
when the daily energy produced by the PV and WT power genera- 
tion plants could only partially cover users’ energy demand. 

A more accurate simulation model of the system is used to 
verify the performance of the control system. In these simulation 
models, implemented in Matlab-Simulink and presented in a 
previous work [26], the main operating parameters, such as cur- 
rent, voltage, pressure and temperature, are evaluated as the 
dynamic behavior of each device with a time step of 1 s. Moreover, 
the verification phase of the EMS was carried out by using actual 
values for solar radiation and wind speed instead of the corre- 
sponding weather forecast. 

Fig. 7 shows the power produced by the two RES systems 
evaluated by using actual data of solar radiation and wind speed. 
The electrical energy produced during the summer week is 


340.0kWh whereas the overall energy requirement is 
268.8 kW h in all cases. A high energy surplus also occurs during 
the spring week (325.0 kWh of overall energy production). The 
weekly energy production by RES diminishes in autumn 
(267.5 kWh) when the difference between energy production 
and demand is very low. A high energy deficit occurs during the 
winter week when the weekly energy produced is 243.1 kWh. 
Fig. 8 shows the charge status of the two storage systems during 
the periods studied, expressed in terms of SOC of batteries and 
hydrogen storage level (the latter is the ratio between the current 
and the maximum hydrogen tank pressure). Since the nominal 
power output of the PV generator is significantly higher than that 
of the WT generator, the energy production of the power genera- 
tion section is greatly affected by the season. In particular, PV pro- 
duction frequently exceeds the energy required by the load during 
the summer so that large amounts of electrical energy are available 
for storage. In this situation, high battery SOC values are often 
reached and therefore the EMS activates the electrolyzer, which 
operates near its nominal condition. Therefore, the hydrogen 
produced is stored and subsequently used to produce electricity 
during periods of great energy deficit as in winter. Unfortunately, 
the current hydrogen storage capacity of the microgrid is not 
enough for long-term storage. 

A better match between energy production and energy demand 
occurs during the autumn and spring weeks. In this case, batteries 
are preferred by the proposed EMS to balance energy production 
and consumption due to their lower utilization costs compared 
to the hydrogen storage system. 

During the winter, electricity production by RES is usually very 
low owing to low solar irradiation and operation of the fuel cell is 
preferred to the batteries to satisfy the large energy deficit in the 
middle of the day. In this case, the fuel cell can work near its 
nominal point with high efficiency whereas a decrease in battery 
efficiency would occur. Therefore, a lower fuel cell utilization cost 
is obtained while a high battery utilization cost (which linearly 
depends on the power supplied, as shown by Eq. (13)), is avoided. 
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Fig. 6. Generation scheduling of hydrogen storage system for the summer case (a), autumn case (b), winter case (c), spring case (d). 
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Fig. 7. Power generation and load demand for the summer case (a), autumn case (b), winter case (c) and spring case (d). 
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Fig. 8. Time evolution of battery SOC and H; tank pressure. 


5. EMS comparative analysis (1) Perfect EMS is based on optimal generation scheduling pro- 
duced by a perfect forecast of wind speed, solar radiation 

To highlight the capabilities of the proposed EMS, in this section and energy demand. In this study, perfect EMS is calculated 
its operational schedules are compared with those of two different by using real data as input values with only one scenario. 


EMS: the so-called “perfect” EMS and the “SOC-based” EMS. Perfect EMS obviously achieves the best results in terms of 
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Fig. 9. Weekly energy flows in the two storage system devices. 


both minimization of utilization cost and maximization of 
overall system efficiency. In real applications, load and 
meteorological data cannot be predicted exactly and perfect 
EMS is used only as a reference and ideal EMS. 


(2) SOC-based EMS is based on optimal generation scheduling 


produced if the utilization costs of batteries, electrolyzer 
and fuel cell are neglected. This EMS considers only battery 
SOC and hydrogen storage level to determine the power 
flows between the different devices. In particular, when gen- 
erated power exceeds the power demand and SOC is below 
100%, the batteries are charged. When the batteries are fully 
charged and the excess power is above the electrolyzer min- 
imum operating power, it is activated until the maximum 
pressure inside the hydrogen tanks is achieved. Vice versa, 
when generated power is less than the power required by 
the user, the batteries start to discharge until a minimum 
value of SOC (SOCyin) is reached. If the batteries reach the 
minimum SOC and there is enough hydrogen in the tanks, 
the fuel cell is turned on. The SOC-based EMS is very simple 
and suitable for real-time control since it requires very low 
computational times. However, this EMS is unable to ensure 
optimal management of the system because it does not take 
into account costs associated with the on/off status of the 
devices. Furthermore, the difference in efficiency with which 
the devices will operate is neglected. 


160 


Figs. 9 and 10 and Table 4 summarize the main results in terms 
of operating hours and utilization costs of the comparative analysis 
of the three aforementioned EMS. 

Results of the proposed EMS should improve with a decrease in 
rolling time and more frequent meteorological data updating 
approaching those obtained by the perfect EMS. However, a 
decrease in the maximum allowable computational time may 
occur and the EMS may not find the optimal solution in time for 
every rolling time. 

On the other hand, the proposed EMS demonstrates several 
advantages in comparison to the SOC-based EMS, especially in 
the summer, winter and spring weeks. Although the total energy 
exchanged with the hydrogen storage system is almost the same 
(Fig. 9), a decrease in utilization cost of the hydrogen storage sys- 
tem occurs in the proposed EMS (Fig. 10). Effectively, with this 
EMS, the hydrogen storage system operates closer to its nominal 
conditions and therefore closer to its maximum efficiency condi- 
tions. In fact, as highlighted in Table 4, an increase in energy con- 
version efficiency is recorded in the proposed EMS with respect to 
the SOC-based EMS even though the results obtained with the per- 
fect EMS could not be reached. In addition, an average decrease in 
the operating cost of about 15% occurs for the proposed EMS com- 
pared to the SOC-based EMS. This reduction in operating cost is 
mainly due to the decrease in the operating hours of the hydrogen 
storage system. 
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Fig. 10. Weekly utilization costs of storage systems. 
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Table 4 
Summary of results obtained from comparative analysis. 


SOCrn (%) Pu, rin (bar) Average ny, (%) Average "parr (%) Load shedding (kW h) Operating cost (€/week) 

Summer 

Proposed 92.6 13.22 51.2 76.2 - 114.4 
Perfect 93.5 13.52 51.4 76.5 - 117.8 
SOC-based 83.9 11.75 38.3 75.6 - 138.4 
Autumn 

Proposed 65.1 8.03 34.9 77.5 - 66.9 
Perfect 82.0 10.32 41.2 77.5 - 61.4 
SOC-based 84.4 10.06 32.2 77.4 - 69.8 
Winter 

Proposed 75.7 4.69 40.4 77.3 31.4 86.5 
Perfect 78.2 3.54 46.3 77.3 27.5 78.1 
SOC-based 72.1 2.86 35.8 76.8 26.7 116.1 
Spring 

Proposed 84.8 10.85 42.7 87.5 - 96.1 
Perfect 80.5 11.20 46.9 87.6 - 94.2 
SOC-based 84.8 10.61 37.8 87.4 - 105.8 


Finally, Table 4 shows that final values of battery SOC and 
hydrogen tank pressures are quite similar for the compared 
EMS. At the end of the summer week, the hydrogen tank pres- 
sure increases by more than 3 bar with respect to the initial 
value (10 bar). This obviously leads to a different initial situation 
for the following week when the storage tank may be 
completely filled. Vice versa, the hydrogen tank pressure reaches 
a lower level at the end of the winter week. This can lead to a 
significant increase in undelivered power for the following few 
weeks. Nevertheless, we wish to emphasize that the purpose of 
this work is to show the benefits and advantages of the pro- 
posed EMS instead of analyzing and optimizing the storage 
capacity of the microgrid studied (the latter will be the task of 
a following paper). 

The main disadvantage of the proposed EMS is the high compu- 
tational time required to find the optimal solution. The use of lin- 
ear approximations to avoid nonlinear equations has led to the use 
of efficient algorithms for linear problems and the certainty of 
finding the global optimum. On the other hand, it increased the 
number of constraints and variables and the complexity of the 
problem. For this case study the problem was solved within 1 h 
on a 2.57 GHz Core i7 processor. The computation time decreases 
significantly to about 20 min for the perfect EMS and 5 min for 
the SOC-based EMS. In all cases, the optimum is found within the 
available decision time (1 h). Further, more accurate solutions with 
a better match with the results of the perfect EMS could be 
achieved by decreasing the rolling time and the update time of 
the wheatear forecast. However, both operations require major 
computational resources. 


6. Conclusion 


In this work, a novel energy management system is proposed 
for a stand-alone microgrid integrated with an energy storage sec- 
tion based on a battery bank and a hydrogen storage system. Start- 
ing from forecasts of weather conditions and load requirements, 
the proposed EMS calculates the optimal generation scheduling 
to minimize operating costs and maximize system efficiency. A sta- 
tistical approach based on the generation of a scenario tree was 
introduced to account for forecasting uncertainties. To compare 
the results obtained with the standard EMS widely applied in the 
literature, a comparative analysis was carried out. In particular, 
the comparative analysis was carried out with reference to four 
different climatic situations (a summer week, an autumn week, a 
winter week and a spring week). 


The results show an improvement in system performance 
compared to the EMS based only on the use of the SOC as a control 
variable. Furthermore, performance deviated only slightly from 
that obtained with a perfect EMS, where input data are perfectly 
forecast. The proposed EMS allows a reduction in operating costs 
of the energy storage section, especially during the summer and 
winter weeks. The increase in the time horizon to more than a 
few weeks and especially the analysis of the annual performance 
of the microgrid can provide useful information on the optimiza- 
tion of storage capacity. 
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